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ABSTRACT 

Accurate and precise photometric redshifts (photo-z’s) of Type la supernovae (SNe la) can 
enable the use of SNe la, measured only with photometry, to probe cosmology. This dra¬ 
matically increases the science return of supernova surveys planned for the Large Synop¬ 
tic Survey Telescope (LSST). In this paper we describe a significantly improved version 
of the simple analytic photo-z estimator proposed by Wang (2007) and further developed 
by Wang, Narayan, and Wood-Vasey (2007). We apply it to 55,422 simulated SNe la gen¬ 
erated using the SNANA package with the LSST filters. We find that the estimated er¬ 
rors on the photo-z’s, + ^phot), can be used as filters to produce a set of photo- 

z’s that have high precision, accuracy, and purity. Using SN la colors as well as SN la 
peak magnitude in the i band, we obtain a set of photo-z’s with 2 percent accuracy (with 
0-(2phot - 2spec)/(l + ^spec) = 0.02), a bias in 2;phot (the mean of Zphot - Zspec) of -9 x 10“®, 
and an outlier fraction (with | (zphot — -Zspec) /(I + -z^spec) I > 0.1) of 0.23 percent, with the re¬ 
quirement that + -^phot) < 0.01. Using the SN la colors only, we obtain a set of 

photo-z’s with similar quality by requiring that + 2phot) < 0.007; this leads to a set 

of photo-z’s with 2 percent accuracy, a bias in Zphot of 5.9 x 10“''^, and an outlier fraction of 
0.32 percent. 

Key words: cosmology; observations, distance scale 


1 INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmological standard 
candles is a corner stone of modem cosmolo gy, and led to the 
observational discover y of cosmic acceleration dRiess et alJI 19981 ; 
IPerlmutter et alj [1999(1 . The true nature of cosmic acceleration, 
dubbed “dark energy” for convenience, remains shrouded in mys¬ 
tery. We do not even know if it is an unknown form of energy 
(hence appropriately named “dark energy”), or a consequence of 
the modification of Einstein’s theory of general relativity (“modi¬ 
fied gravity”)[3 Probing the true nature of cosmic acceleration re¬ 
quires the measurement of both cosm ic expansion history and the 
growth history of large scale str ucture ( iKnox. Song. & Tvsonl2006l ; 
IGuzzo et "nil^OOSl ; IWandlTOO^ . The use of SNe la remains one of 
the most important methods for measuring cosmic expansion his¬ 
tory. 

We can expect a dramatic increase in the number of SNe la 


* E-mail: wang@ipac.caltech.edu _ _ 

^ For recent revie ws, _ see _ iRatra & Vogele-^ j200^; 

FriemaniJ[urnei y & HutereJ |200^: ICaldwel^^^amionkows^ 2005); 
UzaiTl20lBhlwmdd2010ll ; lLi et all 1201 ll) ; IWeinbergerain201^ . 


that can be used to measure cosmic expansion history in the com¬ 
ing years and decades. Thous ands of SNe la are ex pected from the 
Dark Energy Survey (DES) iBernstein et al.|[^012l) . and hundreds 
of thousands of S Ne la are expected from Large Synoptic Survey 
Telescope (LSST) lAbell et alj200^ PI 

The majority of the SNe la from LSST will only have pho¬ 
tometry, and will not have spectroscopic redshifts, due to the ex¬ 
pensive resources required to obtain follow-up spectroscopy. For¬ 
tunately, the redshifts of the SNe la can be estimated using multi¬ 
band photometry (these approximate redshifts are called photomet¬ 
ric redshifts, or photo-z’s). However, the observed photometric SNe 
need to be classified first using photo metry (see “Supernova Pho¬ 
tometric Classification Challenge” by iKessler et al.l J2010bh for a 
detailed discussion). Even applying an optimized version of the top 
performing method from the Supernova Photometric Cl assification 
Challe nge, the photometric-classification algorithm of ISako et al.l 
( I 2 OI ih result in over 25% o f the resultant photom etric SN la sam¬ 
ple remaining non-la SNe dCampbell et alJ[2013l) FI We will defer 


^ http://www.lsst.org/lsst/science _ 

® The photometric-classification m ethod oflCampbell et alj 12013h is based 
on the SN classification technique of lSako etalTuoTlh . but made additional 
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the difficult task of photometric classificationo future work, and fo¬ 
cus on the relatively easier task of estimating redshifts of known 
SNe la using photometry only. 

SNe la with sufficiently accurate and precise photo-z’s can be 
used for probing cosmology. This would dramatically increase the 
science return of supernova surveys planned for LSST. Accurate 
photo-z’s can also enhance the ability of observers to accurately 
target high redshift SNe la for spectroscopy, allowing flexibility 
and optimizatio n in survey design fo r Stage III and Stage IV dark 
energy projects (I Albrecht et al.ll200^ . _ 

In this paper, we build on the pr e vious work by I Wand l l2007ti 
and [Wang, Naravan, & Wood-Vasevi l l2007ti . and present a sim¬ 
ple analytic photo-z estimator. We apply this photo-z estimator 
to 55,422 simulated SNe la generated using the SNANA package 
dKessler et al.l20093l with the LSST filters. We present our method 
in Sec[2 describe the simulation of LSST SNe la in Sec[3 show our 
results in Sec|4] and discuss our findings in Sec|3 


2 THE METHOD 

The a nalytic photo -2 estimator for SNe la proposed by IWand 
(l2007h is empirical, model independent (no templates used), and 
uses observables that reflect the properties of SNe la as calibrated 
standard candles. It was develo ped using the SN l a data released by 
the Supernova Legacy Survey jAstier et al.ll2006h . 

This method requires a training set of SNe la with spectro¬ 
scopic redshifts. The SNe la in the training set are the same as the 
data, but with known redshifts. The training set is randomly chosen 
from the data to span the ranges of the colors and magnitudes of the 
entire data set. This mimics what should be done in practice when 
real data become available: choose a training set from the data to 
carry out spectroscopic redshift measurements. 


2.1 The Prototype Photo-z Estimator 

The prototype photo-z estimator we use here was developed in 
Wang (2007) and Wang, Narayan, and Wood-Vasey (2007). This 
estimator uses the fluxes in griz (or riz) at the epoch of i maxi¬ 
mum flux to make an effective K-correction to the i flux. The first 
estimate of redshift is given by 

2phot = Cl -I- C2gf + C3rf + C4,if + csZf + Cei} + crif (1) 

where gf — 2.51og(/t,), r/ = 2.51og(/r-), if = 2.51og(/i), and 
Zf — 2.51og(/z), and fg, fr, fi, fz are fluxes in counts, normal¬ 
ized to some fiducial zeropoint, in griz at the epoch of i maximum 
flux. 

Next, the photo-z estimator calibrates each SN la in its esti¬ 
mated rest-frame using 

Ails =2.5 log(/fV/i), (2) 

where is the i band flux at 15 days after the i flux maximum 
in the estimated rest-frame, corresponding to the epoch of = 

15(1 -I- 2phot) days after the epoch of i flux maximum. 

The final estimate for the photometric redshift is given by 


-2phot — 


-E 


CiPi, 


(3) 


cuts aided by host-galaxy redshifts with 0.05 < z < 0.55. This enabled 
them to achieve the SN la classification efficiency of 70.8%, with only 3.9% 
contamination from core-collapse SNe. 


where the data vector p = {1, gf ,rf ,if, Zf,rj ,i^, Aii^}. The co¬ 
efficients a (i=l,2,...,8) are found by using a training set of SNe la 
with griz (or riz, for which cg = 0) lightcurv es and measu red 
spectroscopic redshifts. The jackknife technique (lLuDtonlll993h is 
used to estimate the mean and the covariance matrix of a. 


2.2 Photo-z Estimator Applied to LSST Filters 

For application to SNe la from LSST with photometry only, we 
have simplified and modified the analytic photo-z estimator for eas¬ 
ier application and better performance as follows: 

(1) Instead of using the fluxes in all filters except i at the epoch of 
i maximum flux to make an effective K-correction to the i flux, we 
just use the peak magnitudes in all filters for each SN to construct 
an indicator of the spectral shape of the SN. 

(2) We add 2nd order terms in the colors for improved accuracy and 
precision of the estimator. 

(3) We replace Airs with Aii 2 = 2.5 log(//^‘*//i), to increase 

the sample size, where where is the i band flux at 12 days af¬ 
ter the i flux maximum in the estimated rest-frame, corresponding 
to the epoch of = 12(1 -|- Zp^ot) after ths epoch of i 

flux maximum. _ 

Note that our use of Airs was first introduced bv lWangl l l2007h 
to reduce the scatter in the estimated SN la photometric redshifts; 
it is analogous to the use of Amis in SN la calibration, but the use 
of i band m agnitude allo ws us to include a larger set of SNe la in 
our analysis l lWandl2()0% . 

The key to obtaining a set of reliable photo-z’s for SNe is to 
exclude the SNe that are inferred to have unreliable photo-z esti¬ 
mates based on the photometry data alone. We make a first cut in 
SNe to be used as follows: 

(1) Exclude the few SNe that have no i-band photometry. This is 
necessary since we use the i magnitude at maximum light in esti¬ 
mating the photo-z’s. 

(2) Exclude SNe that have measured flux uncertainties at 
maximum light in any passband that exceeds 0.2 magni¬ 
tude. This is useful in reducing the number of outliers with 

K^phot 2:spec)/(l “1“ Zspec )| > 0.1. 

(3) Exclude SNe that do not have i-band lightcurve data that ex¬ 
tends to 12 days after maximum light in the estimated SN rest- 
frame. This enables the refinement of the photo-z estimated using 
the fluxes at maximum light alone. 

We then apply the photo-z estimator as follows: 

(1) Divide the SNe according to the passbands in which the SN 
magnitudes at maximum light are available. 

(2) Eor each set of SNe with photometry in the same passbands, 
construct the initial photo-z estimate given by 

n — 1 

^phot = Cl + ^ ^ Cj.f 1 (rrtp^j TTipj-fi') 

1=1 

n—1 

+ ^ ^ Cn+j(mp^j — rrip^jj^i) 

1=1 

-|-C2n*p -f C2n-|-l*p^ + C2n+2ip^ (4) 

where rrip^j is the magnitude at maximum light in the j-th passband, 
ip is the magnitude at maximum light in the i band, n is the number 
of filters or passbands with available photometry, and {cj}, j = 
1, 2,..., 2n -f 2, are the coefficients to be estimated from using the 
training set of SNe with spectroscopic redshifts. 

(3) Next, calibrate each SN la in its estimated rest-frame using 


(c) 0000 RAS, MNRAS 000, 000-000 


























Photometric Redshift Estimator for SNe la from LSST 3 


Aii2 = ii2d — ip (5) 

where ii 2 d is the i band magnitude at 12 days after the i maxi¬ 
mum light in the estimated rest-frame, corresponding to the epoch 
of = 12(1-1- 2pi,ot) days after the epoch of i maximuml 

light. 

The final estimate for the photometric redshift is given by 
2n+3 

.2phot = ^ ^ Cj Pj , (6) 

j = l 

where the data vector p = {1, (mp,i — mp, 2 ), (trip .2 — 

mp,3),(mp,„_i - mp.„),(mp,i - mp,2)^, (mp,2 - 

mp, 3 )^,(mp,„_i - mp,„)^,ip,ip,i®, Aii 2 }. The coeffi¬ 
cients Cj (j = 1, 2,2n - 1 - 3) are found by using a training set 
of SNe la with lightcurves in n passbands and measured spectro¬ 
scopic redshifts. We use the same modified jackknife technique as 
in Wang, Narayan, and Wood-Vasey (2007) to estimate the mean 
and the covariance matrix of cj (see Sec. 0. 

The improvements that we have made to the photo-z 
estimator in this work reduces the estimated errors on the 
photo-z’s, /(I + Zphot), from 2.5% to 2%, compared to 

Wang, Narayan, & Wood-VasevI ( l2007h . More importantly, it re¬ 
duces the bias in Zphot (the mean of Zphot — Zspec), from —1.5 x 
10~® to —9 X 10“® (a reduction factor of 16.67). 


2.3 Photo-z Estimator Using Colors Only 

Since our ultimate goal is to use SNe la with photometry to probe 
cosmology, it is desirable to have a photo-z estimator that does not 
use the i band flux information, to reduce the likelihood of possible 
double counting of distance information. Note that SN la colors are 
also used in fitting for cosmology, one may argue that including the 
i band flux information is no more double-counting than including 
SN la colors. However, the i band flux of a SN la is closely cor¬ 
related with its distance from us, since we are using the brightness 
of SNe la as distance indicators. SN color is also used in the dis¬ 
tance estimation, but only for the correction to the distance, so it is 
a much smaller effect compared to the SN brightness. Finally, the 
estimation of photometric redshifts are made possible by the use of 
colors, thus they are essential ingredients in any photometric red¬ 
shift estimator. Therefore, having an estimator for SN la photomet¬ 
ric redshifts that uses colors only (the minimal ingredients) enables 
a more robust approach to probing cosmology using photometric 
SNe la. 

We arrive at a photo-z estimator using colors only simply by 
omitting the last 4 terms in Eq. Thus the simplified photo-z es¬ 
timator that uses colors only is given by 

n—1 

^phot ~ Cl + ^ ^ Cj+i(rn-p.j nipj-pi) 

1=1 
n — 1 

+ c„+jimp,j - rripj+if (7) 

1=1 

where nipj is the magnitude at maximum light in the j-th passband, 
n is the number of filters or passbands with available photometry, 
and {cj}, j = 1, 2,..., 2n — 1, are the coefficients to be estimated 
from using the training set of SNe with spectroscopic redshifts. 


3 SIMULATION OF DATA 

The number and location of the LSST deep drilling fields are still 
to be determined; for a general discussion, see the LSST Science 
BookQ Assuming a plausible scenario of 10 fields covering 90 
square degrees, we obtained a set of 55,422 simulated SNe with 
photometric data includin g magnitude uncertai nties, generated us¬ 
ing the SNANA package dKessler et alj|2009^ with the LSST fil¬ 
ters. We have assumed 10 years of observations on th e LSST deep- 
drillin g fields. We have used measured rates of SNla dPildav et al.l 
l2008h as input to the simulations, and imposed the requirement that 
three LSST filters have a SNR greater than 5 (which gives two color 
measurements essential for photo-z estimates). 

We have applied 0.005 mag of smearing to the bandpass ze- 
ropoints; this simulates the variations due to the atmosphere. We 
do not include LSST Photon Simulator effects or systematic band¬ 
pass zeropoint shifts at present^ These effects are expected to be 
smaller than the variations due to the atmosphere. 

To exclude incomplete or very noisy data, we make a few ba¬ 
sic cuts for data quality. Only 55,414 of the 55,422 simulated SNe 
have peak i mags, and only 49,993 have i band lightcurves. Requir¬ 
ing that the i lightcurve data include info of the i peak and extend 
to 12 days in the estimated SN restframe, and that all the peak mag¬ 
nitudes have errors less than 0.2 magnitudes, we arrive at a sample 
of 29,702 simulated SNe. The data quality cuts that we have made 
here are less strict that those made by current surveys; nevertheless 
we find them to be sufficient because of the higher quality photom¬ 
etry expected from the LSST. We will use this sample of 29,702 
SNe la that have passed our data quality cut in the remainder of 
this paper. 

We now describe our Type la SN light curve simulations 
in greater technical detail. We employ the SNANA package 
jKessler et al.ll20099i to simulate the light curves using the SALT2 
Type la SN model. 0 The SNANA package was fi rst used by 
the L SST collaboration to forecast SN observations l l Abell et^ 
12009), and has been us ed extensively by the SDSS collabora¬ 
tion I Kessler et alj2009an. to foreca st observations for the Dark En¬ 
ergy Survey''jBemstd^2^[20l3), and for many other systematic 
studies for supernova cosmology. The specific S ALT2 model used 
is the extended version of the Guy 2010 model i IGuv et al.ll2007L 
bold) . The SALT2 model rest-frame flux F(p, X) as a function of 
phase p and wavelength A, is given in equation[8l The two spectral 
time series, Mq and Mi, and the color law CL are part of the spe¬ 
cific SALT2 input model. The parameters a:o, a:i, c, are randomly 
drawn from distributions that have been fit previously to Type la SN 
data including the effects of dust extinction. Given these parame¬ 
ters the rest-frame flux as a function of phase is available for use. 
The redshift is likewise randomly drawn from a previously mea¬ 
sured distribution, and the observer-frame flux is integrated with 
the LSST filter transmissions. 

F(p,A) =2:o[Mo(p,A)-Fa;iMi(p,A)le-"^-^^^>- (8) 

Our sample of simulated LSST SNe la is representative of data 
expected from the LSST, based on our current knowledge of SNe 


http://www.lsst.org/lsst/scibook 

^ The LSST Photon Simulator (see http://lsstdesc.org/WorkingGroups/PhoSim) 
includes a very detailed simulation of the LSST hardware, and takes a 
star or galaxy location and brightness and propagates the light through the 
optics and camera. 

® The SALT2 model we are using is the most recently updated model avail¬ 
able, with the latest data. 
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la. The SNANA simulations of SNe la have been heavily tested 
on SPSS photometric data up to 2 ; < 0.7 (see, e.g.. lHlozek et aP 
( l2012h : ICamDbell et alj l [20l3 ) , and on SNLS data up to 2 < 1.0. 
This is similar to the redshift range we are using. 


4 RESULTS 

4.1 Photo-z’s using SN la colors and i band peak flux 


We use a slig htly modihed version of the jackknife technique 
lLupton|[l993ll to estimate the mean and the covariance matrix of 
Cj. From the training set containing N SNe la, we extract N sub¬ 
samples each containing N — 1 SNe la by omitting one SN la. The 
coefficients (j — 1, 2,..., 2n -|- 3) for the s-th subsample are 
found by a maximum likelihood analysis matching the predictions 
of Eq. ll6ll with the spectroscopic redshifts. 

The mean of the coefficients Cj (j = 1,2,..., 2n-|-3) are given 
by 




,(«: 

'3 


(9) 


Note that this is related to the usual “bias-corrected jackknife esti¬ 
mate” for Ci, Ci , as follows: 

c/ = cf + (W - 1) (cf - (c,)) , (10) 

where are estimated from the entire traini ng set (with N 
SNe la). [Wang. Naravan. & Wood-VasevI l l2007h found that for 
small training sets (with N < 20) that include SNe la at 2 ~ 0, cj 
give biased estimates of Cj by giving too much weight to the SN la 
with the smallest redshift. For training sets not including nearby 
SNe la, (cj) and cj are approximately equal. We have chosen to 
use {cj) from Eq. as the mean estimates for c, to avoid biased 
2Dhot for SNe la at 2 close to zero JWang. Naravan. & Wood-VasevI 
1200 ^ . 

The covariance matrix of Cj (j = 1, 2,..., 2n -|- 3) are given 
by 


N 

Cov(Ci, Cj) = ^ - (Ci)^ - (Cj)^ (11) 

S = 1 

Using 2phot = CjPj, we find that A2phot = JjpjAcj, since 
the uncertainty in 2phot is dominated by the uncertainty in Cj. 
Therefore estimated error in 2phot is 

{ 2n-|-3 2n-|-3 'j 

EE [pi]Cov{ci,Cj)[pj]\ , ( 12 ) 

i=i j=i J 


where pi and pj (i,j = 1,2, ...2n -|- 3) are the i-th and j- 
th components of the data vector where the data vector p = 

{1, (rrip,! -mp,2), (mp,2-mp,3),..., (mp,„_i -mp,„), (mp,i - 
trtp,2) , (trip,2 tripes) ,.. ., irfip^n — i TTip^n) , ip, ip, ij , At 12 } . 

The simulated SNe are divided into groups according to their 
photometric properties, see Table [T] Each training set for a group 
of SNe consists of the first 100 SNe in that group, assuming that 
the simulated data are randomly ordered. For groups that contain 
very small numbers of SNe, we use the first 20 SNe in each group 
as the training set, or increase the number of SNe used in the train¬ 
ing set until no negative 2phot results. This led us to a training set 
containing a total of 1040 SNe used for fitting the analytic photo-z 
estimator, and a “test set” of 28,662 SNe that are not used in the 


fitting, and are used as a validation set for the analytic photo-z esti¬ 
mator. 

Table[T]shows that we need 1000 SN la spectra for the train¬ 
ing sets for the 29,702 simulated SNe la that satisfy our minimal 
data quality cuts (see Sec.3). When the real LSST photometric SN 
la data are available, we can choose the training set in each color 
group described in Table[T]in a random fashion, while ensuring that 
the SNe la chosen for the training set spans the required color and 
magnitude range. Then spectroscopic redshifts can be obtained for 
the SNe la in the training sets. 

Even 1000 SN la spectra will require substantial observational 
resources. The sizes of the training sets can be reduced to meet 
the constraints of the available resources, which may degrade the 
accuracy and precision of the photometric redshifts. Flowever, the 
performance of the photometric redshifts can be recovered when 
the training sets are expanded later via SN la host galaxy redshift 
measurements. Two strategies can be used to reduce the number 
of SN la spectra needed for the training sets: (1) Remove the SNe 
la belonging to the color groups containing very small fractions of 
the overall sample (G16, G22, G61, and perhaps G21 in Table[T](. 
Note that these minority color groups require a much larger fraction 
of spectroscopic redshifts, and have less well estimated photo-z’s. 
Dropping these will save resources (100 less spectroscopic red¬ 
shifts) with minimal impact on the cosmological constraints. (2) 
Reduce the sizes of the training sets of the other color groups as 
needed. 

In applying the photo -2 estimator to simulated data, we find 
that for each group of SNe with a given set of available photomet¬ 
ric passbands, applying a single fitting formulae to its training set 
may lead to substructure in the estimated 2phot. Therefore, we have 
divided each group into subgroups as shown in Table[T] and derived 
the photo-z estimator for each subgroup using its own training set. 

Table [T] tabulates the resultant precision and accuracy of our 
photo-z estimator for each color group. We measure precision with 
^■(^phot — Zspec)/(1 + 2spec, and the accuracy by the mean of 
2phot — Zspec, as well as the fraction of outliers (also listed are 
the number of outliers that are flagged as having large uncertainties 

itl 2phot)- 

Eig[T]shows (2phot -2:spec)/(l-l-2:spec) versus 2spec for the an¬ 
alytic photo-z estimator applied to simulated SNe from the LSST, 
with cut + 2phot) < 0.1. The cut is made for better 

presentation; only 12 SNe out of 1040 of the training set, and 109 
out of 28662 of the test set have ct/{1 + 2phot) ^ 0.1. Dif¬ 
ferent colors indicate SNe with different photometric passbands, 
as described in Table [T] For SNe represented by the same color, 
different point types represent subdivisions based on brightness to 
improve precision. The histogram of the number of SNe versus 
(2phot — 2spec)/(l + 2spec) Corresponding to Fig[T]is shown in 

FiglSl 

Fig[3]and Fig|4]are the same as Fig[T]and Fig[2 but for impos¬ 
ing the cut -|- 2phot) < 0.01. There are 45 outliers (with 

K^phot - 2spec)/(l -h «spec)| > 0.1) in the bottom panel in Fig[3] 
(the majority are from the rizY color groups Gll through G15), 
representing 0.23% of all 19,640 SNe la from the test sample that 
passed this cut. 

Table [2] shows that we can arrive at photo-z’s with higher 
accuracy, precision, and purity by using the estimated errors on 
2phot to exclude SNe with large + •^phot). If we ex¬ 
clude SNe with -f -Zphot) ^ 0.01, we arrive at a pre¬ 

cision of CT(2phot — 2spec)/(l + 2:spec) = 0.02, a bias in 2phot (the 
mean of 2phot — 2spec) of —9 x 10“®, and an outlier fraction (with 
|(2phot - 2spec)/(l -I- 2spec)| > 0.1) of 0.23 percent (see FigO. 
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Figure 1. The performance of the analytic photo-z estimator applied to sim¬ 
ulated SNe from the LSST, as indicated by (Zphot ~ ^spec)/(l + ^spec) 
versus Zgpec, tvith cut + 2phot) < 0.1 for cleaner presenta¬ 

tion (the numbers of outliers are given in Table[T).. Different colors indicate 
SNe with different photometric passbands, as described in Table[T] For SNe 
represented by the same color, different point types represent subdivisions 
based on brightness to improve precision. 


Note that for the threshold of -|- Zphot) < 0.01, there 

are no outliers in the training set. We can require that there are no 
outliers in the training set in order to arrive at this threshold for 

'^Zphot/(l + 2phot). 


4.2 Photo-z’s using SN la colors only 

As we discussed in Sec l2.3l having an alternative photometric red- 
shift estimator that uses SN la colors only is helps increase the 
robustness of cosmological constraints obtained using photomet¬ 
ric SNe la. Table[T]tabulates the resultant precision and accuracy of 
our photo-z estimator using SN la colors only, for each color group. 
The colors only results are tabulated immediately below the results 
for colors plus i band peak flux (see the previous subsection), and 
indicated by “Y” under the column for “colors only”. Using SN 
colors only leads to slightly worse precision, but similar accuracy 
for the resultant photo-z’s. 

FigOand Fig|7]are results corresponding to Fig[T]and Fig[3 
but for Zphot obtained using Eq.l]?), which uses SN la colors 
only. Using the SN la colors only, we obtain a set of photo- 
z’s with similar quality to using SN la color and i band peak 
flux by requiring that + .^phot) < 0.007; this leads 

to a set of photo-z’s with 2 percent accuracy, a bias in Zphot of 
5.9 X 10“"^, and an outlier fraction of 0.32%. The threshold level 
of + Zphot) < 0.007 is set naturally by requiring that 
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Figure 2. The distribution of (Zphot — Zspec)/(1 + Zspec) from Fig[T] 


there are no outliers in the training set. There are 58 outliers (with 
|(Zphot - Zspec)/(1 + Zspec)| > 0.1) in the second panel from 
the top in Fig|7](the majority are from the rizY color groups Gll 
through G15 — similar to Fig[^, representing 0.32% of all 18,370 
SNe la from the test sample that passed this cut. We do not expect 
such a small fraction of outliers to have noticeable impact on the 
cosmological fit; we will investigate this in future work. 


5 SUMMARY AND DISCUSSION 


We have shown that a simple analytic photo-z estimator for SNe 
la can be built for the LSST, given suitably chosen training sets of 
SNe la with spectroscopic redshifts that are representative of the 
general properties of SNe la fr om the LSST . While our method 
is based on that proposed by ( l2007h . and advanced by 

[Wang. Naravan. & Wood-VasevI ( 2007 ), it represents a significant 
improvement over previous work (which require the use of both 
colors and i magnitudes of SNe la) in enabling the use of SN la col¬ 
ors only. This eliminates the possible complications that can arise 
from using SN la brightness in estimating photo-z’s, since the peak 
brightness of a SN la is closely correlated with its distance from us. 

Our sample of simulated LSST SNe la is representative of data 
expected from the LSST, based on our current knowledge of SNe 
la. The multiple-band photometry from the LSST enables us to di¬ 
vide the SNe la into different sets (see Sec l2.2l and Table [T}. For 
each set, a photo-z estimator is constructed using only quadratic 
functions of colors (see Eq.(7)), with the constant coefficients given 
by fitting to a training set that contains a maximum of 100 SNe la 
with spectroscopic redshifts (see Table [T](. Note that the division 
of SNe la into color groups is empirical, dependent on the proper- 
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Figure 3. The same as Fig[3 but with cut cTzphot /(^ ^phot) < 0-01 for 
improved accuracy, precision, and outlier reduction. 


ties of the training set, assumed to be representative of the overall 
sample. 

In addition to how the photometric SNe la should be grouped, 
and the performance of our photo-z estimator for each color group, 
Table^contains other helpful information as well. As expected, the 
performance of the photo-z estimator improves with the increasing 
number of filters in which photometric data are available; note that 
knowing the filters in which the SNe la are visible provides a crude 
estimate of the redshift range of the SNe la. Not surprisingly, the 
performance of the photo-z estimator generally degrades with the 
increasing redshift of the SNe la. 

While Table^lists all the SNe la that passed our data quality 
cut (see Sec[3}, Table[^shows us how to derive a high purity sample 
of photo-z’s by making cuts on the SNe la based on the estimated 
errors of the photo-z’s. This is very useful since we are systematics 
limited, rather than statistics limited, in using photometric SNe la 
for cosmology. 

The performance of our simple analytic SN la photo-z es¬ 
timator is quite impressive when applied to simulated LSST SN 
la photometric data (see Table [Hi. We find that the estimated er¬ 
rors on the photo-z’s, -1- aphot), can be used as filters 

to produce a set of photo-z’s that have high precision, accuracy, 
and purity. Using SN la colors as well as SN la peak magnitude 
in the i band, we obtain a set of photo-z’s with 2 percent accuracy 
(with a{zpi,ot - 2spec)/(l -b ^spec) = 0.02), a bias in ^phot (the 
mean of Zphot — «spec) of —9 x 10“®, and an outlier fraction (with 
|(Zphot - Zspec)/(1 -b Zspec)| > 0.1) of 0.23 percent, with the re- 




Figure 4. The distribution of (^phot — ^spec)/(1 -b Zspec) from Figjs] 
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Figure 5. The same as Fig^ but for Zp^iot obtained using SN la colors 
only. 
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Figure 6. The distribution of (2:phot — ^spec)/(l + 2spec) from Fig|^ 


quirement that + -^phot) < 0.01. Using the SN la colors 

only, we obtain a set of photo-z’s with similar quality by requiring 
that + Zphot) < 0.007; this leads to a set of photo-z’s 

with 2 percent accuracy, a bias in 2phot of 5.9 x 10“"^, and an outlier 
fraction of 0.32 percent. 

In future work, we will investigate whether the precision and 
accuracy achieved by our photo-z estimator is adequate for cosmol¬ 
ogy, We will first need to simulatte more realistic samples of LSST 
SNe, which include non-la’s. We will need to expand and advance 
our photo-z estimator for weeding out non-la SNe. There is much 
work to do before our photo-z estimator can be applied to the real 
LSST data. 

In principle, our method for estimating photo-z’s should be ap¬ 
plicable to other surveys, but its performance depends on the preci¬ 
sion of the photometry of the surveys. However, it can still be used 
as a quick way of estimating photo-z’s that is model-independent, 
and complementary to the template-based methods. The applica¬ 
tion of our method to non-la SNe and galaxies can also be ex¬ 
plored. This method for estimating photo-z’s for SNe la was ini¬ 
tially inspire d by a similar method for estim ating galaxy redshifts 
presented bv lWang. BahcalL & Tume3(ll998l) . It would be interest¬ 
ing to develop a simple photo-z estimator for LSST galaxies based 
on this method. 

We expect that our SN photo-z estimator for the LSST will 
be useful in probing dark energy using LSST SNe la with pho¬ 
tometry only. It will also provide an independent cross-check 
for SN photo-z’s d erived from other, template-based photo-z es ¬ 
timators (see, e.g., iKirn & Miauell (l2007h ; iKessler et all ( l2010ah : 
IPalangue-Delabrouillel ( |201(]|) 1. Our linear photo-z model is very 
convenient for propagating photometric uncertainties, especially 
compared to more “black-box” type machine learning algorithms. 
In future work, we will examine the cosmological constraints that 
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Figure 7. The same as Fig|^ but with cut -F ZpUot) < 0.007 

for improved accuracy, precision, and outlier reduction. 


can be obtained using LSST SNe la with photometry only, with 
photo-z’s estimated using the method presented in this paper. 


ACKNOWLEDGMENTS 

We are grateful to Michel Wood-Vasey and Alex Kim for helpful 
discussions, and Kirk Gilmore for an internal review of our paper 
on behalf of the LSST Dark Energy Science Collaboration. YW 
was supported in part by NASA grant 12-EUCLID12-0004. 


REEERENCES 

Abell, P. A. et al. 2009, e-prints arXiv:0912.0201 
Albrecht, A., et al.. Report of the Dark Energy Task Force; available 
online at arXiv:astro-ph/0609591vl 
Astier, R, et al. 2006, A & A 447, 31-48 (2006) 

Bernstein, J.P. et al., Astrophysical Journal 753 (2012) 152 
Caldwell, R. R., & Kamionkowski, M., 2009, 

Ann.Rev.Nucl.Part.Sci., 59, 397 
Campbell, H., et al., 2013, ApJ, 763, 88 
Dilday, B., et al., 2008, ApJ, 682, 262 
Frieman, J., Turner, M., Huterer, D., ARAA, 46, 385 (2008) 
Guzzo, L., et al., 2008, Nature, 451, 541 
Guy, J. et al. 2007, Astronomy & Astrophysics, 466, 11 
Guy, J. et al. 2010, Astronomy & Astrophysics, 523, A7 
Hlozek, R., et al., 2012, ApJ, 752, 79 


© 0000 RAS, MNRAS 000, 000-000 



























































8 


Wang, Gjergo, & Kuhlmann 
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Table 1. Division of SNe into sets according to available photometric pass- 
bands, and into subsets in each set for improved precision. The “outliers” 
column indicates the number of SNe in the training set and test set that have 
I (^phot — ^spec)/(l + -Sspec)] > 0.1. The “flagged in test” column indi¬ 
cates the number of outliers in the test set that have cTzphot /(^ ^phot) ^ 
0.1. 
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Table 2. Obtaining sets of photo-z’s with higher accuracy, precision, and 
purity by using the estimated errors on 2phot to exclude SNe with large 

^-Sphot /(I + ^phot)- 
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